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ABSTRACT 

Angular differential imaging is a high-contrast imaging technique that reduces 
quasi-static speckle noise and facilitates the detection of nearby companions. A 
sequence of images is acquired with an altitude/azimuth telescope while the in¬ 
strument field derotator is switched off. This keeps the instrument and telescope 
optics aligned and allows the field of view to rotate with respect to the instrument. 

For each image, a reference PSF is constructed from other appropriately-selected 
images of the same sequence and subtracted to remove quasi-static PSF struc¬ 
ture. All residual images are then rotated to align the field and are combined. 
Observed performances are reported for Gemini North data. It is shown that 
quasi-static PSF noise can be reduced by a factor ~5 for each image subtraction. 

The combination of all residuals then provides an additional gain of the order of 
the square root of the total number of acquired images. A total speckle noise at¬ 
tenuation of 20-50 is obtained for one-hour long observing sequences compared to 
a single 30s exposure. A PSF noise attenuation of 100 was achieved for two-hour 
long sequences of images of Vega, reaching a 5-sigma contrast of 20 magnitudes 
for separations greater than 8". For a 30-minute long sequence, ADI achieves 
30 times better signal-to-noise than a classical observation technique. The ADI 
technique can be used with currently available instruments to search for ~ 1 
Mj up exoplanets with orbits of radii between 50 and 300 AU around nearby 
young stars. The possibility of combining the technique with other high-contrast 
imaging methods is briefly discussed. 
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1. Introduction 

Direct detections of very faint exoplanets and brown dwarfs near bright stars are essential 
to understand substellar formation and evolution around stars. This endeavor is now one 
of the major goals for next generation 10-m telescope instruments and future 30- to 100- 
m telescopes. The task is dauntingly difficult. The exoplanet or brown dwarf image is 
usually much fainter than the background from the brilliant stellar point spread function 
(PSF) image. Besides the Poisson noise, ground-based telescopes suffer from atmospheric 
turbulence that produces random short-lived speckles that mask faint companions. If these 
two limitations were the only ones, a simple solution would be to integrate longer to average 
these random noises and gain as the square root of the integration time. But observations 
have shown that, for integrations longer than a few minutes, the PSF noise converges to a 
quasi-static noise pattern, thus preventing a gain with increasing integration time (Marois et 
al. 2003, 2005; Masciadri et al. 2005). To achieve better detection limits, it is thus necessary 
to subtract the quasi-static noise using a reference PSF. Both ground- and space-based 
imaging are plagued with this stellar PSF calibration problem caused by imperfect optics 
and slowly evolving optical alignments. For ground-based imaging, subtraction of a reference 
PSF obtained from a star close to the target achieves a factor ~ 4 of PSF noise attenuation, 
leaving residuals that are also quasi-static and thus severely limiting detection of fainter 
companions (Marois et al. 2005). For space telescopes that have a better PSF stability, like 
HST, a partial solution was found by subtracting two stellar images acquired during the same 
orbit with a different roll angle. This technique, called “roll deconvolution”, successfully 
subtracts the stellar image by a factor 50 (Schneider & Silverstone 2003; Fraquclli et al 
2004) but is also ultimately limited by PSF evolution. A similar technique, called angular 
differential imaging (ADI), can be used on ground-based altitude/azimuth telescopes to 
subtract a significant fraction of the stellar quasi-static noise and can potentially achieve 
detection limits that improve as the square root of the integration time. 


1 Based on observations obtained at the Gemini Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the 
Gemini partnership: the National Science Foundation (United States), the Particle Physics and Astronomy 
Research Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), the 
Australian Research Council (Australia), CNPq (Brazil) and CONICET (Argentina). 
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In this paper, the ADI technique is described and its performance is analyzed using a 
simple analytical model and data from the Gemini Altair adaptive optics system and the 
NIRI near infrared camera. The PSF stability with Altair/NIRI is studied and its impact 
on ADI performances is discussed. Detection limits for three stars of our ongoing young 
nearby star survey are then shown. A comparison between ADI and normal imaging is also 
presented. Finally, the possibility of using ADI with other high-contrast imaging techniques 
is discussed. 


2. The Angular Differential Imaging Technique 

ADI is a PSF calibration technique that can, in principle, suppress the PSF quasi¬ 
static structure by a large factor (Marois 2004). It consists of the acquisition of a sequence 
of images with an altitude/azimuth telescope and the instrument rotator turned off (at the 
Cassegrain focus) or adjusted (Nasmyth) to keep the instrument and telescope optics aligned. 
This setup improves the stability of the quasi-static PSF structure throughout the sequence 
while it causes a slow rotation of the field of view (FOV) with respect to the instrument. 
Note that only the FOV, not the PSF, rotates with time. Since the FOV rotates during 
an exposure, companion PSFs are smeared azimuthally. This effect increases linearly with 
angular separation and slightly decreases companions peak intensity. Short exposures and 
the use of an optimized aperture photometry box can minimize this effect. For each image, 
after data reduction and image registration of the whole sequence, a reference PSF obtained 
from other images of the same sequence is subtracted to remove the quasi-static structure. 
Given enough FOV rotation during the sequence, this subtraction preserves the signal from 
any eventual companion. All the image differences are then rotated to align the FOV and 
are median combined. 

This technique offers a number of advantages over more classical ground-based obser¬ 
vations since the target observations themselves are used to construct a reference PSF. This 
means that the reference PSF has the same spectrum and brightness as the target and that 
no time is lost to acquire reference observations of a different target. Ghost images from 
optical reflections and the sky flux are also removed by the subtraction. The detector flat 
field errors are averaged since the FOV is integrated with different pixels as it rotates on the 
detector. The ADI technique is a generalization of the “roll deconvolution” technique used 
with HST since several images, each at a different field angle, are acquired and combined. 
A technique similar to ADI has been developed independently by Liu (2004) to search for 
circumstellar disks. 
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3. Noise Attenuation Theory with ADI 

For each reduced and registered image of an ADI sequence, a reference PSF has to be 
built from the images of the same sequence. The way that this reference PSF is built is of 
great importance since it directly affects the noise attenuation performance. We have used 
two methods to construct the reference PSF. 

The first method is simply to take the median of all the images of the sequence. If 
enough field rotation has occurred during the sequence so that an eventual point source has 
moved by at least twice its full width at half maximum (FWHM), then this point source 
will be largely rejected by the median which will leave only the average PSF. The minimum 
radial separation at which this occurs is noted R mm . Since the median is taken over a large 
number of images, the pixel-to-pixcl noise (i.e. PSF, flat field, dark and sky Poisson noises 
and detector readout noise) of the reference image is much less than that of any individual 
image. Thus this first method minimizes the noise in regions where the residuals are limited 
by pixel-to-pixcl noise. However, since a sequence typically lasts more than an hour, the 
reference PSF only has modest quasi-static speckle correlation with the individual images of 
the sequence. 

The second method is to take the median of only a few images as close in time as 
possible but for which the displacement due to field rotation between the images is, at a 
given separation, at least 1.5 PSF FWHM. This displacement ensures that the flux inside 
the PSF core of an eventual point source is not significantly reduced by the subtraction. The 
time r min required for such field rotation is function of the separation angle from the target, 
the target azimuth A and zenith distance z and the telescope latitude qb. The rotation rate 
ifj (degree/minute) of the FOV is obtained from the time derivative of the parallactic angle 
and is given (McLean 1997) by 




0.2506 COS ' 4W 
sin ^ 


( 1 ) 


Figure 1 is provided as a reference to determine the time interval for observations from the 
summit of the Mauna Kea. This second technique provides better quasi-static speckle noise 
attenuation since the reference PSF is built using images acquired at short time intervals. 
However, the pixel-to-pixcl noise of the reference image may not be negligible compared to 
that of an individual image. 

The first method is optimized for regions where the residuals are limited by pixel-to- 
pixcl noise while the second is optimized for regions where the residuals are still limited by 
speckle noise. The combination of both techniques into a single reduction algorithm will be 
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discussed in section 5.2. 

The ADI technique attenuates the PSF noise in two steps: (i) by subtraction of a 
reference image to remove correlated speckles and (ii) by the combination of all residual 
images after FOV alignment to average the residual noise. The noise attenuation N/AN is 
defined as the ratio of the local noise N in an image over the noise AN of the residual image. 

The noise attenuation obtained by the subtraction of the reference PSF, [A/AA] S (#, r, f exp ), 
is a function of the angular separation, 9 , the time interval r between the image and its refer¬ 
ence and the individual image exposure time f exp (including overheads). Strong quasi-static 
speckle correlation between successive images leads to strong attenuation and thus better 
detection limits for a given total integration time. A good and stable seeing is thus expected 
to deliver better ADI quasi-static speckle attenuation. For a perfect case when all static 
and quasi-static speckles have been removed by ADI, detection limits are ultimately limited 
by short-lived atmospheric speckles that have a correlation timescale of a few tens of ms, 
shorter than the time interval required to obtain sufficient FOV rotation to build the ADI 
reference PSF. 

The additional noise attenuation resulting from the combination of n de-rotated residual 
images is function of the correlation of those images. This attenuation is given by ^/n//g, 
where n e g is the effective number of uncorrelated residual images; if all residuals are uncor¬ 
related, then n e g = n. Equivalently, for an image sequence of total time t = t exp n, one has 
to wait for a time 


n 

Tdcorr 

^eff 


exp 


( 2 ) 


before two de-rotated residual images are decorrelated. Thus the total noise attenuation 
N/AN will be 


N 
A N 


(f, 9 , T, f e x P ) 


rsj 


Tdcorr 


N 
A N 


(9, t, t, 


exp; 


JS 


( 3 ) 


For simplicity, let us consider a dominant speckle noise having a correlation timescale 
T sp eck- The ADI noise attenuation resulting from the combination of all difference images is 
function of T spe ck- Its behavior can be defined for two limiting regimes: a) when either r spe ck 
is much longer than r or shorter than t exp and b) r spec k is longer than t exp but shorter than 
r. 


In the first regime, the residuals of consecutive image differences are decorrelated, either 
because the correlated structure of the PSF which lasts for long times (r spec k t) has been 
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removed, leaving only Poisson noise (rd CO rr = t exp ), or because the PSF structure was already 
uncorrelated between consecutive images in the first place (T spec k < t exp =>■ r dcorr = t exp ). 
In this case, the speckle attenuation increases with the square root of the total number of 
image differences. The decorrelation timescale is simply 


^dcorr ^exp* (4) 

In the second regime, since T spec k is longer than t exp but shorter than r, the residuals of 
consecutive differences are partially correlated and do not average-out as efficiently. Assum¬ 
ing that only the noise having a spatial scale ~ 1 PSF FWHM limits point source detection, 
the time needed to decorrelate this noise is the shortest time between t F whm, the time for a 
one FWHM FOV rotation, and r spec k. The decorrelation timescale is thus 

Flcorr — MIN(t F WHM, Apeck)- (5) 


We note that if r is several minutes and if individual ADI differences are dominated 
by random short-lived speckles, then the noise attenuation can follow Eq. 4 for several 
images but ultimately converges to Eq. 5 when sufficient averaging of short-lived speckles 
has occurred. 

It is important to known that reducing t exp to increase the number of images for a given 
image sequence of length t would not result in better detection limits. For the first regime, 
shorter t exp would mean a smaller companion S/N per image; the final S/N would thus be 
the same (neglecting losses due to overheads). For the second regime, a shorter f exp would 
yeild quasi-static speckles correlated over more consecutive images, also resulting in the same 
performance. 

We emphasize that ADI guarantees a gain in detection with increasing observing time 
for both regimes. In the worst case, the speckle attenuation efficiency is limited by field 
rotation. This is an enormous advantage over classical observation in which the quasi-static 
aberrations prevent significant gain after a relatively short observing time. An equation 
similar to Eq. 5 also applies to classical observations, but in this case there is no tfwhm 
since the field is not rotating. Speckle attenuation is then governed by r spec k, which can 
be very large, reducing drastically the efficiency of the observations. Section 6.4 presents a 
comparison of ADI and classical observations. 
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4. Observations 

The ADI technique was first used at the Gemini North telescope using the Altair adap¬ 
tive optic system (Herriot et al. 1998) and the near-infrared camera NIRI (Hodapp et al. 
2000) in queue mode. Observations and detailed results obtained for three stars of our nearby 
young star survey, Vega, HD18803 and HD97334B are presented. Data for two other stars, 
HIP18859 and HD1405 will also be discussed for comparison since they have been acquired 
with a different technique. The exposure time was fixed to 30s, a time short enough to 
minimize PSF saturation and the off-axis PSF smearing clue to the FOV rotation during the 
exposure and long enough to get a good observing efficiency and to be limited by the sky 
background rather than by the read noise at wide separations. 

For Vega, HD18803 and HD97334B, the imaging sequence consisted in the acquisition 
of a series of saturated images in the off methane 1.58 /jrn 6.5% bandwidth filter with the 
Cassegrain rotator turned off. This filter was chosen to minimize the brightness ratio of the 
star to that of methanated exoplanets or brown dwarfs. Unsaturated short exposures were 
acquired before and/or after the saturated sequence to calibrate detection limits. 

For Vega, data were obtained on 2004 August 26 and 2004 September 1 (program 
GN-2004A-Q-11). In total, 225 and 177 30s exposures were respectively acquired on the 
26th and September 1st. Because of Vega’s brightness, short exposure PSFs of a nearby 
reference star were acquired for photometric calibration. Seeing conditions were excellent 
on the 26th and average on the 1st. Images were saturated inside a ~6" diameter. For 
HD18803, 90 30s exposures were obtained on 2004 December 24 (program GN-2004B-Q4) 
and, for HD97334B, 90 30s exposures were acquired on 2005 April 18 (program GN-2005A- 
Q16). For HD18803 and HD97334B, seeing conditions were average to good and PSFs were 
saturated inside diameters of 0.7". Table 1 summarizes the observations. Strehl ratios were 
obtained by analyzing unsaturated data acquired before and after each sequence and by 
comparing the PSF peak intensity with that of a simulated unaberrated PSF having the 
same pixel sampling, bandpass and integrated flux. It can be deduced from table 1 that 
since R min is less than the saturated radius for all three targets, ADI can be applied at all 
separations to detect point sources. 

As part of Altair science verification we have obtained a sequence of observations of the 
star HIP18859 on 2003 November 18 (program GN-2003B-SV-102). During this sequence, 
the filter was switched from the broadband H filter to a narrow band filter every fourth 
exposure to acquire unsaturated images. This data set will be discussed in section 6.1. 

Observations of the star HD1405 were obtained on 2004 August 23 (program GN-2004B- 
Q-14) with the instrument rotator operating to keep the FOV orientation fixed throughout 



the sequence. In total, 38 30s exposures were obtained. These observations will be used in 
section 6.4 for a comparison between ADI and classical observations. 


5. ADI Data Reduction Algorithm 
5.1. Preliminary Data Reduction 

The data reduction consists of flat held normalization, bad pixel correction using a 
median over surrounding pixels, and distortion correction using software provided by the 
Gemini Staff (Trujillo, private communication) and modified to use the IDL interpolate 
function with cubic interpolation. Images were then copied into larger blank images to 
ensure that no FOV was lost when shifting and rotating images. The center of the PSF of 
the first image of the sequence was then registered to the image center by minimizing the 
diffraction spikes residuals after subtraction of a 180-degree rotation of the image. The rest 
of the images were then registered by cross-correlation of the diffraction spikes with the first 
image. An azimuthally symmetric intensity profile was finally subtracted from each image 
to remove the smooth seeing halo. 


5.2. ADI algorithm 

As discussed in section 3, two methods can be used to subtract the quasi-static PSF 
structure: subtracting the median of all images or subtracting a reference PSF obtained from 
a few images acquired as close in time as possible. These two methods can be combined into 
a single algorithm that optimizes speckle subtraction and minimizes pixel-to-pixel noise. 
First, the median of all the images is subtracted from each individual image. For a sequence 

h(tl, $l)> -^2(^2 5 $2 )i -^3(^3) $ 3 )) • • • j 7 n (f n , @n) (6) 

of n reduced and registered images (see section 5.1), where is the mean time of exposure 
i and is the FOV orientation at time tj, the first reference subtraction is simply 

/f = Ii- median(I : , I 2 ,1 3 ,..., I n ) (7) 

An optimized reference PSF (second method) is then obtained for each image by median 
combining 4 images (two acquired before and two after) that show at least a 1.5 FWHM (one 
FWHM in the off methane 6.5% filter is equal to ~ 3 pixels) FOV orientation difference. 
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This choice insures that the average r of the reference PSF is ~0 and that the linear time- 
evolution of quasi-static speckles having T spec k > 2r is removed. For the construction of this 
reference PSF, the image is broken into many annuli to accommodate for the dependence of 
7 ~min on the separation. The intensity of the reference PSF is then scaled appropriately inside 
each annulus to minimize the noise after subtraction. The scaling factor converges to zero if 
the annulus is dominated by pixel-to-pixcl noise or to unity if it is dominated by correlated 
speckles. The optimized reference PSF is then subtracted. This step can be summarized by 
the following equation for each annulus 

l t DI = l ? - « x median(l[i b , 1° b _n 1° c , I? + c+1 ) (8) 

where a is a normalizing factor to minimize the noise inside the annulus and b and c are the 
number of images required to get at least a 1.5 FHWM FOV orientation difference between 
It and the images acquired respectively before and after Jj. Differences are then rotated to 
align the FOV to that of the first image. Finally, a median is taken over all differences. The 
final combination step is thus 


Ip DI = median [lf DI , rot(I^ DI , A6 >i_ 2 ), rot(I^ DI , A6 »i_ 3 ), ..., rot(I^ DI , A6>i_ n )] (9) 

Optionally, the final residual image may be convolved by a Gaussian of FWHM equal to that 
of the PSF to attenuate further the high frequency noise. Table 2 summarizes the entire 
ADI reduction algorithm. 


6. Results 

6.1. PSF Evolution Time-Scale 

The PSF noise evolution timescale can be studied through the evolution of the noise 
attenuation [iV/AA] for the difference of two images as a function of the time interval, r. 
The noise attenuation is the ratio of N and AN, which are the measured rms noises in 
an annulus of width equal to one PSF FWHM in the original and in the difference images 
respectively. For this analysis, all images have first been unsharp masked using an 8 x 8 
FWHM box to remove low-frequency noise and then median filtered with a 1 x 1 FWHM 
box to remove hot/bad pixels. This step is necessary to prevent biasing the noise estimate 
N of single images and leaves only speckles that have a spatial scale of the order of one 
FWHM. Images were subtracted two by two with increasing time interval. Fig. 2 shows the 
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noise attenuation for Vega, HD18803 and HD97334B for angular separations of 2, 4, 6 and 
8 ". 


All noise attenuation curves show stronger noise attenuation for shorter time intervals. 
For HD18803 and HD97334B, the noise attenuation reaches ~ 2 — 5 inside 4" for r ~2 min¬ 
utes. At larger separation, the residuals are limited by pixcl-to-pixcl noise. For separations 
less than 4", the noise attenuation drops by a factor 2 after approximately 60 minutes. For 
Vega, the noise attenuation reaches ~ 4 — 6 at all separations for r ~ 1 minutes and drops 
by a factor 2 after approximately 15 to 20 minutes. The stronger noise attenuation achieved 
on Vega and HD97334B for short time intervals could be explained by better seeing that 
stabilizes the structure and enables a better subtraction. Additionally, since r min is less than 
10 minutes at V within ±lh from the meridian for all targets having a declination —30° 
to +65° at Mauna Kea (see Fig. 1), ADI can be used with little loss of speckle attenuation 
(~30%) on most of the available sky for separations greater than ~ 0.5". 

Analysis of the HIP18859 data set that included a frequent filter change to acquire 
unsaturated images shows that a drop by a factor of 2 in speckle attenuation occurs following 
each filter change. This evolution of the PSF structure is probably clue to the filter wheel 
not returning to its exact position after each change. This suggests that the best observing 
procedure is to prevent any alteration of the optical setup during an observing sequence to 
maximize the PSF noise stability. 


6.2. ADI Speckle Attenuation 

Fig 3 shows for all ADI targets the average noise attenuations [A/AA]g achieved for 
one image difference and the noise attenuation [N/ A N] obtained after median combining all 
the image differences. Again, an 8x8 FWHM unsharp mask and a 1 x 1 FWHM median 
filter were applied to each image to produce this figure. Each ADI difference attenuates by a 
factor 2-10 the quasi-static speckle noise. Note that these attenuations are better than what 
is observed in Fig. 2 since here the optimized reference PSF is the result of a median of 4 
images, 2 acquired before and 2 after. Random speckles are averaged and the linear speckle 
evolution with time of speckles having at least r spe ck A 2r m i n is removed. Similarly to what 
was observed in Fig 2, better seeing conditions (HD97334B compared to HD18803) seem to be 
related to better speckle noise attenuation. A total noise attenuation of ~ 35 between 0.8 and 
10" was obtained for HD18803 and HD97334B, while this attenuation reached 100 for Vega 
(August 26th data). The higher attenuation for Vega comes partly from the larger number 
of images (225 vs 90) and partly from better seeing conditions which provided a better 
attenuation from single image subtractions. HD97334B and HD 18803 noise attenuations are 
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shown to improve at all separations down to the detector saturation limit and clearly below 
it if we extrapolate the performances shown at 0.8". 

The bottom panel of Fig. 3 is generated by dividing the ratio of the noise attenuation 
of the combined ADI difference [Af/AA/] over the noise of a single ADI differences [AT/AAT]g 
by the expected noise attenuation y 7 n if the noise was decorrelated. From Eq. 3 



Because the pixels affected by the secondary mirror support diffraction spikes are masked 
from the images, the effective number of images combined is reduced at some separations. 
The bottom panel of Fig. 3 has been corrected for this effect, which can reach 20% at small 
separations. For separations greater than 2", all three targets achieve more than 70% of the 
a Jn attenuation expected if residuals are decorrelated. At smaller separations the attenuation 
is lower, revealing a correlated noise in successive ADI differences that lowers the number of 
independent images (see Eq. 5). 

Since ADI does not achieves a a Jn noise attenuation gain at all separations, the bottom 
panel of Fig. 3 and Eq. 10 can be used to estimate the residual speckle decorrelation time that 
is currently limiting ADI performances. Fig. 4 upper panel shows the estimated decorrelation 
timescale as a function of angular separation. Typical decorrelation times are of the order 1-3 
minutes and are generally smaller at larger angular separations. Normalizing these curves by 
the time MAX(t FWH m, t exp ) (see Fig. 4 bottom panel), it can be deduced that, for separations 
less than 3", since Td CO rr is less than tfwhm, ADI differences are limited by residual speckles 
having r spe ck < Tfwhm- These residual speckles evolve faster than the time needed to obtain 
an ADI reference PSF and cannot be subtracted by ADI. For HD97334B and HD18803, for 
separations greater than 3", ADI is limited by the FOV rotation or the noise is decorrelated 
(limited by t exp ). For Vega, the decorrelation time increases until tfwhm = 4x P and then 
decreases. The fact that Td CO rr is significantly greater than tfwhm when tfwhm > f exp suggests 
that Vega ADI differences, at those separations, are limited by correlated noises that are 
bigger than one FWHM. This noise thus requires more than a FWHM FOV rotation to be 
decorrelated. The upper limit at ~ 1.5 FWHM is expected since a minimal FOV rotation 
of 1.5 FWHM (r m in) has been chosen to build the ADI optimized reference PSF. The ADI 
optimized reference PSF subtraction is a temporal filter that guarantees that no noise can 
have a decorrelation time longer than the time needed for a ~ 1.5 FWHM FOV rotation. 
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6.3. Contrast Performances 

Fig. 5 shows detection limits (5a) in magnitude difference as a function of angular 
separation obtained with the ADI technique for all three ADI targets. To produce this 
figure, data were reduced following the procedure explained in section 5.1. No unsharp mask 
was used for the ADI reduction since multiple tests have shown that while this filter is 
effective for suppressing the low-frequency spatial noise, it does not improve candidate S/N 
and the photometry is slightly biased in the process. The detection limits are calculated 
using the ratio of simulated companion peak intensities over the noise in the residual image 
as a function of angular separation. The flux normalized unsaturated PSF was then used 
to simulate smeared companion PSFs as a function of angular separation. The noise is 
calculated inside annuli of increasing diameter and width equal to 1 PSF FWHM. To account 
for the PSF smearing effect due to FOV rotation, the image and simulated companions were 
convolved with elliptical Gaussian of one FWHM in the radial direction and one FWHM 
plus a smearing term in the azimuthal direction, the smearing term ranging from zero at the 
center to typically two FWHM at 10", depending on the rotation rate of the FOV. Finally, 
detection limits were corrected for the Altair estimated anisoplanatism following the Strchl 
S equation found in the Gemini web page 4 

S(6) = 5 0 e"(^) 2 (11) 

where 6 is expressed in arcsec. The detection limit obtained with the ADI technique on 
Vega at separations greater than ~5" is two orders of magnitude deeper than the Palomar 
H -band image (Metchev et al. 2003) and approximately a factor of ten deeper than the Keck 
AT-band image of (Macintosh et al. 2003). The current speckle attenuation (20-100) achieved 
with ADI is comparable or better to what is currently obtained (50) from one HST orbit 
using the roll subtraction technique. Although the ADI technique is inherently optimized for 
relatively large separations, the good Gemini PSF stability enables excellent performances 
at sub-arcsec separations. Indeed, the ADI contrast of Am = 11.1 — 11.9 (5<r) at 0.8" 
obtained on HD18803 and HD97334B equals the Am = 11.0 (5<r) at 0.8" obtained with the 
simultaneous differential imaging (SDI) VLT camera optimized for multi-wavelength speckle 
suppression (Biller et al. 2006). This shows the potential of ADI to achieve high-contrast 
detection at sub-arcsec separations using a simple, yet efficient, observing technique with 
standard instruments. Fig. 6 illustrates the noise attenuation obtained for Vega (August 
26th) using the ADI technique. 

4 http://www.gemini.edu/sciops/instruments/altair/altairlndex.html 
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Mass limits corresponding to these observations, corrected for the filter use 5 , are esti¬ 
mated using evolutionary models of Baraffe et al. (2003) assuming ages of 350, 45 and 85 Myr 
for Vega, HD18803 and HD97334B, respectively (Song et ah 2001; Montes et ah 2001). Both 
HD 18803 and HD97334B achieve detection limits of 1-2 Mj up at 3" (60 AU for both targets), 
while ~3 Mj up is obtained for Vega at 8" (63 AU). The ADI technique is thus well suited 
to survey jovian companions at intermediate separations (50-300 AU) orbiting young nearby 
stars. 


6.4. Comparison between ADI and Classical Imaging 

In the previous sections it was shown that the ADI technique can achieve high contrast 
given a sufficiently long integration time and good PSF stability. To compare the perfor¬ 
mances of ADI and classical observations we analyze the first 38 images of the HD97334B 
ADI sequence and the 38 images of the HD1405 “classical” sequence. For this analysis, both 
data sets have been reduced according to section 5.1. Furthermore, an 8 x 8 FWHM unsharp 
mask was applied to all images to remove the low spatial frequency quasi-static noise. Then 
a 1 x 1 FWHM median filter was applied to all images to remove the bad/hot pixels. These 
steps are performed here only to bring the classical observations on even ground with ADI 
in order to study the evolution of the noise at spatial scales that most severely limit point 
source detections. For the HD97334B sequence, images differences were obtained according 
to section 5.2; these differences were then rotated to align the FOV to that of the first image. 

An increasing number of images (differences for HD97334B) of both sequences were 
median combined to study the noise attenuation as a function of the total observing time at 
2", the results are presented in Fig. 7. The ADI reduction technique achieves 30 times better 
speckle noise attenuation compared to classical AO observations in 30 minutes integration 
time. This figure also illustrates the power of ADI imaging in which noise attenuation, and 
thus companion S/N, increases nearly as the expected A Jn while it saturates rapidly for 
normal imaging. 

5 The star-to-planet brigthness ratio is reduced by a factor ~2.6 when using a 6.5% bandpass methane 
filter instead of a broad band H, as derived from a theoretical spectrum of Allard et al. (2001). 



7. Discussion 


ADI is a general high-contrast imaging technique that can be applied to any existing or 
upcoming large altitude/azimuth telescope. It is also flexible enough to be combined with a 
number of other high-contrast imaging techniques. 

ADI at small separations (< 1") requires long time intervals and thus may suffers from 
PSF variations. These variations come from variable seeing and slowly evolving quasi-static 
aberrations from the telescope and instrument optics. Unsaturated data inside 0.8" are 
required to estimate performances. Extrapolation from our observations predicts that noise 
attenuation will be limited by speckles evolving faster than Tfwhm required for a 1 FWHM 
FOV rotation. Speckle noise attenuation of the order of a factor 10 should be feasible for 
one-hour sequences. The use of a multi-wavelength instrument (Marois et al. 2000; Doyon et 
al. 2004; Lafreniere et al. 2004; Marois et al. 2004, 2005; Biller et al. 2004) or an IFU (Sparks 
& Ford 2002) to acquire simultaneous images at multiple wavelengths across the methane 
absorption bandhead at 1.6 /mi through the simultaneous spectral differential imaging (SSDI) 
technique (Racine et al. 1999; Marois et al. 2000; Biller et al. 2004; Marois 2004; Marois et 
al. 2005) could provide good short-lived and common-path speckle attenuations and increase 
detection limits. It has been shown that SSDI instruments are ultimately limited by non- 
connnon path aberrations, which are expected to be stable over long periods of time as 
they come almost entirely from the instrument itself (Marois et al. 2005). Hence, the ADI 
technique nicely complements SSDI since it can be used to subtract the residuals caused by 
the non-common path aberrations. 

Future high-contrast instrumentation for 8-10 m class or larger telescopes based on 
high-order adaptive optics (AO) systems (Macintosh et al 2004; Mouillct et al. 2004) will 
most likely improve the stability of the PSF. Thus, if combined to such instruments, ADI 
could prove even more successful. These new AO systems will also provide the high Strchl 
ratios required to bring coronagraphy at the forefront of high-constrast imaging. However, 
even a very good coronagraph cannot totally suppress the light from uncorrected quasi¬ 
static wavefront errors and some level of quasi-static speckle noise will inevitably be present 
in coronagraphic observations. ADI would be a nice addition to coronagraphy as it could 
attenuate those residual speckles. 

Ideally, all the techniques mentioned above, ADI, SSDI, high-order AO and coronag¬ 
raphy, could be used together to form an extremely powerful tool to detect exoplanets and 
brown dwarfs around stars. 



8. Conclusion 


The ADI observing technique was described and its performance using Altair/NIRI at 
Gemini was presented. It was shown that faint companions can be detected with better S/N 
when compared to classical observing techniques for a wide range of declinations. The ADI 
technique produces a reference PSF from the same target imaging sequence, removing the 
need to move to a nearby star for PSF calibration or to acquire sky exposures (for H -band 
imaging). Since the reference PSF is built using images acquired minutes apart, the reference 
PSF shows a good quasi-static speckle correlation. 

The stability of the PSF plays a crucial role in ADI as it not only determines the 
speckle attenuation from the reference image subtraction but it also determines the regime 
in which the noise is attenuated with increasing observing time. It was reported that at 
Gemini with Altair/NIRI using 30s exposures, the PSF evolves on timescales of ~ 10 — 
60 minutes and the attenuation by subtraction of a reference image reaches ~ 2 — 6 for 
short time intervals, achieving better speckle attenuation with better seeing conditions. The 
observations of HIP18859, for which a filter change during the sequence reduced significantly 
the speckle attenuation, underscore the importance of maintaining the optical setup fixed 
during the sequence. It was shown that the gain in S/N with increasing total observing time 
for separation greater than 2" reaches more than 70% of the optimal case, indicating that the 
noise is mostly decorrelated between residual images for these separations. Typical residual 
speckle decorrelation time is of the order of a few minutes. The speckle noise residuals 
decorrelate faster for object having faster FOV rotation. In all cases, ADI guarantees a 
larger gain with longer observation sequences. To our knowledge, this is the first time that 
such behavior is clearly demonstrated for an acquisition and reduction technique designed for 
speckle attenuation. The wall raised by quasi-static speckles that prevents a gain with longer 
integration time for standard observing techniques (Marois et al. 2003, 2005; Masciadri et 
al. 2005) can thus be removed by ADI. Comparison with a classical imaging technique shows 
that ADI achieves 30 times better speckle attenuation in 30 minutes integration time. 

The noise attenuation obtained on Vega was 100, reaching a contrast of ~20 magnitudes 
at 8" separation (63 AU). Observations of the young stars HD 18803 and HD97334B yielded 
detection limits in difference of magnitude of 11.1-11.9 at 0.8", similar to the SDI camera 
at VLT (Am of ~11 at 0.8"), which is an optimized speckle suppression instrument. When 
combined to substellar models and estimated age for these stars, these observations show 
that ADI is well suited to search for jovian companions having a mass greater than 1-2 Mj up 
50-300 AU away from nearby young stars. Finally, ADI could easily and advantageously be 
combined with SSDI, high-order AO and coronagraphy to improve the detection limits of 
exoplanets and brown dwarfs at all separations. 
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Tabic 1: Observations 


Object 

Date 

Nb 

images 

Field rotation 
(degree) 

-Rmin (") 

r min (min) 
at meridian & 1" 

Strehl 

Vega 

08/26/04 

225 

99 

0.15 

8 

0.24 a 

Vega 

09/01/04 

177 

69 

0.1 

8 

0.12 

HD18803 

12/24/04 

90 

99 

0.15 

2 

0.10 

HD97334B 

04/18/05 

90 

54 

0.3 

6 

0.16 


a Short exposures for the August 26th run are saturated. Since the September 1st short exposures fill 55% 
of the pixel electron well at the PSF peak intensity, it is approximated that the Strehl ratio for August 26th 
is at least 2 times higher than the one estimated for September 1st. 



Table 2: ADI Data Reduction Algorithm 


Reduction 

Flat field normalization 

Bad pixel correction 

Distortion correction 

Copy in bigger blank image 


Processing 

Image registration 

Radial profile subtraction 

Diffraction spikes & saturation masking 


ADI 

Subtraction of the median of all images 
For each 0.6" wide annulus 

Find 4 images with r > 1.5r min 
Median combine these 4 images 
Flux normalization 

Subtract reference annulus 


Calculate parallactic angle 

Rotate images 


Coaddition 

Median combination of all differences 
(optional) Convolution by Gaussian 
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ig. 1.— Interval of time required for a point source to move by 1 X/D (1.6 //rn on a 
arneter telescope) at 1" as a function of hour angle for various declinations. Calculate* 
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Fig. 2.— Noise attenuation obtained by subtracting images two by two with increasing time 
interval for HD97334B, HD18803 and Vega acquired on the 08/26/04 and on the 09/01/04. 
Solid, dashed, dot-dashed and triple dot-dashed lines are for 2", 4", 6" and 8" respectively. 
For Vega, there is no solid line for 2" since images are saturated at that separation. Small 
lines show the estimated noise attenuation limit imposed by photon, sky, flat, read and dark 
noises. 
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Fig. 3.— Single ADI difference [A/AA],, and total ADI [A/A A] noise attenuations with sep¬ 
aration for HD97334B, HD18803 and Vega. Dotted lines show attenuation from a single ADI 
difference while solid lines show the attenuation after median combining all ADI differences. 
The bottom panel shows the ratio of the speckle noise attenuation of the combined ADI 
difference [A/AA] over that of a single ADI difference [A/AA]s normalized by the square 
root of the total number n of images in each sequence. Data for HD18803 are truncated at 
9" due to a PSF decenter that brought a part of the field > 9" outside the FOV. 
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Fig. 4.— Upper panel: Estimated residual speckle decorrelation timescale Td CO rr as a function 
of angular separation for HD97334B, HD18803 and Vega (August 26th). Bottom panel: 
Normalized residual speckle decorrelation time as a function of angular separation for the 
three same targets. The star symbol on each curve indicates at what separation tfwhm = t exp . 
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Fig. 5.— Final ADI detection limits (5<r) as a function of angular separation for HD18803, 
HD97334B and Vega (August 26th). The initial detection limit for the HD97334B PSF 
(dot-dashed line) is also shown to illustrate the ADI noise attenuation performance. 
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Fig. 6.— Vega (August 26th) ADI data reduction. A: A single image after flat field normal¬ 
ization, bad pixel correction, distortion correction, registering and removal of an azimuthally 
symmetric profile. FOV is 22" x 22" using a linear intensity range of ±10 -6 from the esti¬ 
mated PSF peak intensity. B: A single ADI difference image shown with the same intensity 
range. C: same as B with an intensity range 25 times smaller. D: The final combination 
of all ADI differences shown with the same intensity range as C. The central saturated 6" 
diameter region as well as diffraction from the secondary mirror supports have been masked. 
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Fig. 7.— Normalized noise attenuation from the median combination of an increasing number 
of images (differences in the ADI case) at 2 " separation. Each single image represents 30s 
exposure time. The dotted and dashed lines are respectively for the ADI reduction and the 
HD1405 without held rotation sequence. See section 6.4 for more details. 






